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ABSTRACT 

We study the kinetics of two-species annihilation, A + B — > 0, when all particles 
undergo strictly biased motion in the same direction and with an excluded volume 
repulsion between same species particles. It was recently shown that the density 
in this system decays as i -1 / 3 , compared to t -1 / 4 density decay in A+B — ► with 
isotropic diffusion and either with or without the hard-core repulsion. We suggest 
a relatively simple explanation for this t -1 / 3 decay based on the Burgers equation. 
Related properties associated with the asymptotic distribution of reactants can 
also be accounted for within this Burgers equation description. 
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1. INTRODUCTION 



While the kinetics of two species annihilation, A + B — > 0, has been extensively 
investigated and reasonably well understood [1-5], surprising results have been recently 
obtained in one dimension when both species move according to driven diffusive motion 
[6]. In this model, each particle can attempt a move only to the right and the attempt 
is successful only if the particle lands on an unoccupied site. If the particle lands on a 
site already occupied by a member of the same species, the move attempt is rejected and 
the initial particle remains at its starting position. However, if the particle lands on a site 
which is occupied by a member of the opposite species, then an A + B — > reaction occurs 
and both particles disappear. This motion of each species may be viewed as strongly biased 
diffusion with a hard core repulsion between same species particles. 

For this system, one might naively anticipate that, by performing a Galilean transfor- 
mation to a zero momentum reference frame, the long-time kinetics should be the same as 
that of A + B — > with isotropic diffusion. For this latter system, the density is known 
to decay in time as t -1 / 4 when the initial densities of the two species are equal [1-5]. This 
result has been found to hold both for the case of free diffusion and when an excluded- 
volume repulsion between same species particles exists. In fact such a Galilean universality 
has been derived for the case of single species reactions and a special type of two-species 
reaction [6]. 

However, extensive numerical simulations in one dimension indicate, rather surpris- 
ingly, that for two-species annihilation with driven motion, the density asymptotically 
decays as t -1 / 3 for a Poissonian initial distribution of the two species [7]. It was also 
found that the approach of the density to the t -1 / 3 asymptotic form was rather slow. 
Even after 10 6 time steps, local estimates of the density decay exponent were still varying 
systematically with time. On the other hand, for isotropic diffusion, stable estimates of 
the decay exponent are achieved after 10 3 - 10 4 time steps. As we discuss below, there is 
an interesting reason for this slow approach to asymptotic behavior in the driven system. 

In this paper we provide a theoretical argument for the t -1 / 3 decay of the density 
in A + B — > with driven diffusive motion. Our approach is based on first elucidating 
the dynamics of a single AB interface when the two species are initially separated (Sec. 
II). A continuum description in terms of the inviscid Burgers equation, together with 
continuity of particle currents, shows that the velocity of the interface can be determined 
in terms of the relative densities of the two adjoining species. This information is used to 
infer that when the initial condition is homogeneous and random, a typical single-species 
domain grows as t 2 / 3 . This result, together with the fact that the density inside a domain 
is simply related to initial density fluctuations, leads to the t -1 / 3 decay of the density. 
In Sec. Ill, we present numerical simulations of the reaction, focusing, in particular, on 
the density decay exponent, the domain length growth exponent, and the behavior of 
the scaled domain density profile. This latter consideration has previously been found 



2 



to be fruitful in understanding both the kinetics and the spatial structure of reactants in 
A + B — > with unbiased diffusion of each species [8]. Our numerical results support many 
of our qualitative arguments. Finally, in Sec. IV, we provide general conclusions as well as 
discuss the kinetics with arbitrary biases of each species and in higher spatial dimension. 

2. ANALYTICAL APPROACHES 

A. Single species motion 

Let us first recall some basic properties of single-species driven diffusion in one dimen- 
sion where the positions of individual particles are updated sequentially and in random 
order [9]. Denote by rii(t) the occupation number for the i th site at time t. By the 
excluded- volume interaction, rii is restricted to be either or 1. Depending on the local 
configuration of particles on sites i — 1, i, and % + 1, the occupation number for site % at 
time t + 5t has the following three possible outcomes, 



m(t + St) = 



rii(t) with probability 1 — 2aSt, 

ni(t)rii+i(t) with probability aSt, (1) 



, rii(t) + (1 — ni(t))rii-i(t) with probability a8t. 

Here 5t is the microscopic time step and aSt is the probability that a particle attempts to 
move during the time step. These three possibilities have the following origins: If neither 
the i th nor the (i — l) th particles were selected for a move attempt, then rii(t) does not 
evolve, as indicated in the first line. The second line describes the state of site i at t + St 
if a move attempt is performed on the particle initially at site i. Similarly, the third line 
gives the state of site i if a move attempt is performed on a particle initially at site % — 1. 

Corresponding to these elemental updates, the Master equation for the average occu- 
pation number rii(t) is, 

ni(t + 6t) - rii(t) = aSt [rTi(t)n i+1 (t) - ni(t)ni-i(t) + n i - 1 (t) - rii(t)}. (2) 

In writing this equation, we have exploited the fact that for one-dimensional driven dif- 
fusion in a finite ring or in an open system but sufficiently far from the boundaries, the 
occupation numbers for the different sites are uncorrelated [9] . This implies that a mean- 
field description is exact so that the two-particle correlation function, n,i(t)nj(t), may be 
replaced by the product of single-particle densities, ni(t)fij(t), in Eq. (2). In the contin- 
uum limit, the leading order expansion of the above Master equation yields the Burgers 
equation, 

dn Ax d r _._, _. , Ax 2 d 2 n /0 , 
Tt=-MTx Hl - n)]+ 2At^ (3) 
where Ax is the lattice spacing and At = a -1 is a macroscopic time step during which 
one move attempt is performed, on average, for each particle. In the following we choose 
the length and time scales such that Ax = 1 and At = 1. 
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This Burgers equation may be solved by first reducing it to the diffusion equation by a 
Cole-Hopf transformation [10]. Since the latter equation is trivially soluble, the solution to 
the Burgers equation may be obtained by performing an inverse Cole-Hopf transformation 
on the corresponding diffusive solution. A sketch of the Burgers equation solution in the 
long-time limit for a localized initial source n(x, 0) = A5(x) is shown in Fig. 1. There are 
three universal regions in the long-time solution which are independent of the microscopic 
details of the initial density distribution. First, there is a constant-slope "interior" region 
where n(x, t) « (1 — x/t)/2,with x in the domain t — 2\/~At < x < t. This behavior arises 
because the leading low-density particles have a relatively large velocity and continue to 
propagate faster than the denser trailing particles. At either end of this interior region, 
i. e., at x = t — 2\f~At and at x = t, there are two diffusive boundary regions whose widths 
grow as ~ \ft. In the interior region, the solution to Eq. (3) is asymptotically close to 
the profile of the solution to the inviscid Burgers equation (shown dotted in the figure). 
This inviscid solution consists of the straight line segment n\ m \ sc \&(x,i) = (1 — x/t)/2 for 
t — 2\fA± < x < t, while r!i nv i SC id(x, t) = otherwise. 

B. Dynamics of a single interface 

To understand the evolution of domains in the initially homogeneous system it will 
be helpful to determine the dynamics of a single AB interface. Therefore consider the 
"separated" initial condition where n^(x,t = 0) = Ha for x < and n^(x,t = 0) = for 
x > 0, whereas ns(x,t = 0) = for x < and ns(x,t = 0) = for x > (Fig. 2). 
The evolution of such a system provides essential insights into the nature of the reaction 
front [11-16], and has proven useful in understanding spatial structure in A + B — > with 
homogeneous initial conditions, since relatively sharp reaction fronts naturally form as the 
system evolves. For the separated system, we will apply a mean-field approach in which 
the free motion of the A and B particles is modeled by the inviscid Burgers equation. Then 
by invoking mass balance at the AB interface, we determine the interface velocity and the 
spatial distribution of reactants. 

Assume that the interface initially at x = moves from x to x + dx in a time dt. Thus 
the Bs within the slab [x, x + dx] will be completely replaced by As. Consequently, the 
incident flux of A, jA, must be sufficient to annihilate all the Bs which do not escape from 
the slab and also repopulate the slab at the final density of As. That is, 

2a dt = ha dx + R, (4a) 

where R represents the As which will undergo reaction in the slab. Similarly, the outgoing 
flux of B, 2b, includes the initial number of Bs in the slab minus those that disappear by 
reaction: 

2b dt = us dx — R. (46) 
These two lead to the mass balance equation, 

(j A (x,t) +j B (x,t))dt = {n A {x,t) +n B (x,t))dx, (5) 
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where dx = vab dt, with vab being the interface velocity. In the inviscid Burgers equation 
these currents are related to the corresponding particle densities, tia and ng, by, 



3a,b{x,€) = n A ,B(x,t)(l - n A ,B(x,t)). (6) 

The concentration of A at the interface depends upon whether the interface itself 
moves faster or slower than the point x which defines the beginning of the depletion zone 
for the A distribution (Fig. 2). According to the Burgers equation, this boundary point 
moves with velocity v = X — ItLa- Thus if the interface moves to the right faster than x, the 
hypothesized depletion zone exists and the concentration of A at the interface coincides 
with solution of the inviscid Burgers equation. On the other hand, if vab < v, the depletion 
zone does not appear and the concentration of A is constant throughout the domain. Thus 
for the concentration of A at the interface, we have 

r (l-x/t)/2 for vab > v 
n A (x,t) = l (7) 
I riA otherwise. 

On the other hand, since the concentration of B remains constant throughout the domain 

n B (x,t) = n B . (8) 

Substituting Eqs. (6)-(8) and the relations dx/dt = vab and v = 1 — 2ua into the mass 
balance equation (5) gives, 

{1 - 2n s ( v / 2 - 1) for n A >n B (V2-l) 
( 9 ) 
for n A <n B (V2-l). 

To check whether these formulae account for the propagation of an AB interface on a 
microscopic level, we numerically simulated the A + B — > reaction with initially separated 
species which undergo the same driven diffusive motion. The agreement between the 
simulation data and the model predictions indicate that our phenomenological description 
based on mass balance and the inviscid Burgers equation accurately accounts for the 
microscopic behavior (Fig. 3). However, the actual width of the reaction zone is primarily 
controlled by diffusion and therefore its behavior is beyond the scope of an inviscid model. 



C. Density decay and domain growth rates 

Let us now return to the original problem of equal initial concentrations of A and B 
which are randomly distributed. As time increases, a region initially rich in a particular 
species (due to initial density fluctuations) will evolve into a single-species domain. To 
determine the time dependence of the size of such a domain, note that if the domain length 
is L, the initial excess of a particular species and thus the number of particles remaining 
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within it after the minority species is annihilated is typically of order \[L. Consequently, 
the concentration in a typical domain is of order [1-5] 

n(t) ~ l/^W). (10) 

On the other hand, from Eq. (9) the rate of change of the domain length dL/dt is equal 
to the difference of the velocities of its two boundaries, 

dL/dt oc | vab — vba I • (11a) 

Also from Eq. (9) this velocity difference is proportional to the density difference across 
the interfaces. These differences are typically of order the domain densities themselves. 
Thus we conclude that 

dL/dt ~ ua — tib ~ n(t). (11°) 

Combining Eqs. (10) and (11), we predict that the domain size grows as L ~ t 2 / 3 , while 
the density decays as 

n^-t- 1 / 3 . (12) 

It is also instructive to estimate the leading correction to this result. Recall that 
the difference between the solution of the inviscid and viscid Burgers equation occurs in 
the boundary regions whose widths are proportional to ~ \ft (Fig. 1). Consequently, the 
inclusion of a viscosity term should lead to at 1 / 2 correction to the typical domain length, 
i. e., L ~ At 2 / 3 + Bt 1 / 2 . Correspondingly, this gives rise to a correction of order ~ t -1 / 6 to 
the asymptotic behavior of the density. This slowly vanishing correction is a possible reason 
for slow convergence of the previous [7] and our simulations to the predicted asymptotic 
behavior. 



3. SIMULATION RESULTS 

We performed simulations on a system of 5 x 10 5 lattice sites for 1.5 30 = 1.9 x 10 5 time 
steps and averaged over 30 initial configurations. As has been observed previously [7], the 
local estimate of the density decay exponent from our data is systematically decreasing with 
time (Fig. 4(a)). An asymptotic value of —1/3 for this exponent is not incompatible with 
our data. A larger scale simulation would help resolve this issue. Another basic quantity 
which also exhibits slow approach to asymptotic behavior is the average domain length 
L(t) (Fig. 4(b)). The local exponent of this quantity grows systematically with time and 
appears to approach the asymptotic value 2/3 that was predicted in our phenomenological 
approach. 

A deeper understanding of the behavior of the system can be gained by examining the 
distribution of particles within a domain. In the case of A + B — > with isotropic diffusion, 
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we examine both the "micro-canonical" and "canonical" domain profiles [8] . To construct 
the former, we first define the extent of a particular domain of A, say, by the distance 
between the two B particles which enclose this A domain. We then rescale the lengths 
of individual domains to a constant value and then superpose the corresponding rescaled 
densities. Finally, it is convenient to plot the scaled profile p(x/L,t)t 1 ' 3 , where p(x/L,t) 
is the local density within a domain at scaled position x/L and time t (see Fig. 5). These 
data collapse onto a universal curve, except for the leading and trailing regions, where the 
slopes for later time seem to be steepening. We suggest that the departure from scaling can 
be partially understood by recalling that the boundary layers around each domain grow as 
y/i, while a domain by itself grows as t 2 / 3 . Thus the relative width of the boundary layer 
decreases as t -1 / 6 , leading to the slope of the scaled density profile as t 1 / 6 . To illustrate 
this effect, we plot, in Fig. 6, the trailing tail of the microcanonical domain profile when 
the ordinate is rescaled according to t 1 / 6 . The apparent existence of scaling according to 
t 1 / 6 in this subregion suggests that there are really two length scales for a single domain. 
One, growing as t 2 / 3 , corresponds to the "bulk" domain growth and is described by the 
inviscid Burgers equation. The other, corresponding to the boundary layer, grows as ~ \ft 
and may be described only by including the higher-order viscosity term in the Burgers 
equation. 

We have also examined the canonical density profile in which the densities of all do- 
mains are superposed about a common center, but without rescaling the domains lengths 
before the superposition. Strikingly, these profiles are nearly symmetric about a point 
which is centered some small distance upstream from the domain center (Fig. 7) . Further- 
more, the decay of this profile is quite accurately described by a simple exponential. 



To summarize and to place our results in an appropriate context, it is helpful to 
consider the A + B — > reaction when there are arbitrary hopping anisotropies Pa and 
Pb, respectively, for the A and B species. Thus for a given species, the probability for 
a particle to jump to the right is (1 + fi)/2 and to the left is (1 — P)/2. The model we 
have considered thus far corresponds to the extreme case of Pa = Pb = 1, i- e., particles 
hop only to the right. The general case can be usefully described in terms of the "phase 
diagram" whose axes are the degree of anisotropy of each species (Fig. 8). 

In the continuum limit, it straightforwardly follows that for arbitrary but equal biases 
of A and B, (0 < /3a = Pb = P < 1), the Burgers equation, Eq. (3) will be modified only 
by the introduction of the coefficient P which multiplies the first-order spatial derivative, 



4. DISCUSSION 



dn 
~dt 




ld 2 n 



(13) 



2 8x 2 ' 
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Such a coefficient in the equation of motion can be absorbed by a suitable rescaling of x 
and t, leading to universal behavior in the hopping anisotropy, as long as this anisotropy 
is non-zero and the same for both species. Hence the density decay exponent along the 
line Pa = Pb is expected to be equal to —1/3, except for (Pa, Pb) = (0,0). This singular 
point corresponds to the case of pure diffusion where the density decay exponent equals 
— 1/4, independent of the existence of a same species hard core repulsion [1-5,8]. 

The remainder of the phase diagram corresponds to systems with unequal biases for the 
two species. Applying the arguments of Sec. II, the evolution of the typical domain length is 
determined by the difference in the velocities at each interface. This gives dL/dt ~ Pa~Pb, 
so that the average domain size grows linearly in time. Consequently the concentration 
varies as n(t) ~ l/y/L ~ 1/Vi. Our simulations of A + B — > with an unequal biases for 
A and B are consistent with this heuristic prediction. 

Thus diffusion-controlled two-species annihilation with same-species hard-core repul- 
sion has rich variety of kinetic behaviors in one dimension. The density decay exponent 
equals 1/4, 1/3, and 1/2 for zero, equal, and unequal biases, respectively. For the first 
and third situations, rate equation approaches have failed to give the correct asymptotic 
behavior. In contrast, for the equal bias case we have shown that the Burgers equation ap- 
proach, supplemented by mass balance at the interfaces between domains, explains several 
characteristics of the system. More subtle features, such as the density profile inside the 
domains, are beyond the simple-minded approach developed here but might be accessible 
by a refinement of our technique. It is also worth mentioning that macroscopic characteris- 
tics of the system, e. g. : the typical domain size L ~ t 2 / 3 , naturally arise in other contexts 
where the Burgers equation plays the leading role [17]. 

Finally, it is worth considering the effect of driven motion on A + B — > in greater 
than one spatial dimension. Numerically, we have investigated two natural versions: In 
the NE model, we allow particle to either move up or to the right with equal probabilities, 
while in the NES model, we allow particles to move either up, down, or to the right with 
equal probabilities. For both cases, the density appears to decay asymptotically as t -1 / 2 . 
Interestingly, there is a significant time range where the density decay exponent is clearly 
less than —1/2. This same anomaly also appears in A + B — > with isotropic diffusion of 
the reactants, however. This can be ascribed to the existence of a limited temporal regime 
where there is some penetration of particles of one species into a domain of the opposite 
species. This effect is even more strongly pronounced in three dimensions [18]. 

Qualitatively, one can estimate the effect of driven diffusion on the reaction kinetics 
in arbitrary dimension d. If L is the typical length of a single-species domain in the 
drift direction, then its volume V scales as Lt^ d ~ 1 ^ 2 , since the domain length in the 
remaining d — 1 transverse directions still grows by diffusion. Repeating now the steps 
employed in the previous treatment of the one-dimensional situation, we have dL/dt ~ n 
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and n ~ V 1//2 ~ L 1 / 2 t ^ d ^/ 4 . These can be readily solved to yield 

f t -(d+i)/6 for rf < 2 
t -d/4 for 2 < d < 4 



(14) 



for d > 4. 



Our simulations are covered by the first line of Eq. (14). When d > 2, the longitudinal 
length predicted by the above argument, L ~ t^ 5 ~ d ^ 6 , is smaller than the diffusional 
length scale. This indicates that the drift is irrelevant above two dimensions and the 
reaction kinetics is the same as two-species annihilation with isotropic diffusion. 
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Figure Captions 

Figure 1. Asymptotic solution of the Burgers equation for the initial value problem n(x, t = 
0) = AS(x) with ( ) and without (• • •) the viscosity term. 

Figure 2. Schematic illustration of evolution of the interface between two initially separated 
A and B domains. The motion of each species is described by the inviscid Burgers 
equation. There is a depletion zone in the A distribution which begins at x and 
extends to the interface. In a time dt, the interface advances from x to x + dx and 
the initially uniform B distribution in the slab [x,x + dx] is replaced with an slowly 
varying A distribution whose density is approximately equal to the A density at the 
interface. 

Figure 3. Density profile near the AB interface at t = 0(o),t = 50(+),t = 100(*),t = 150(x), 
and t — 200(A), for 300 configurations with ha = ns = 0.5. The calculated positions 
of the interface are marked by the o. The dotted line shows the calculated value of 
riA = 0.207 for the A density at the interface. 

Figure 4. Time dependence of the local exponents for: (a) the density decay exponent, and (b) 
the domain length exponent. 

Figure 5. The scaled micro canonical density profiles at t = 1.5 18 = 1477(A), t = 1.5 21 = 
4987(x),t = 1.5 24 = 16833(*),t = 1.5 27 = 56814(+), and t = 1.5 30 = 191750(o). 

Figure 6. Trailing edge of the microcanonical density profiles when rescaled by the factor t 1 ^ 6 . 
The symbols refer to the same times as Fig. 5. 

Figure 7. The canonical density profile at t = 1.5 21 = 4987. 

Figure 8. "Phase diagram" for the density decay exponent for different hopping anisotropies of 
the A and B species. 
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